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INTRODUCTION 

It has been known for the past four decades that Quantum Chromodynamics (QCD), together with the electroweak 
interactions, underlies all of nuclear physics. Soon after the discovery of QCD it became apparent that the complexities 
of the theory at strong coupling would hinder analytic progress in understanding the properties of the simplest hadrons, 
let alone the simplest features of the nuclear forces. Wilson pointed the way to eventual direct quantitative confirmation 
of the origins of nuclear physics by formulating Lattice QCD [1], a regularization and non-perturbative definition of 
QCD that is suitable for the intensive computational demands of solving QCD in the infrared. Only in the last five years 
or so has Lattice QCD emerged from a long period of research and development — where only qualitative agreement 
between calculations and experiments could be claimed — to the present, where precise predictions for hadronic 
quantities are being made. In particular, presently, fully-dynamical calculations with near-exact chiral symmetry at 
finite lattice-spacing have become standard, with lattice volumes of spatial extent L> 2.5 fm and with lattice spacings 
in the range b% 0.12 fm. It is still the norm that the light-quark masses, m q , are larger than those of nature, with typical 
pion masses m K ~ 300 MeV. 

Nuclear physics is a vast, rich field, whose phenomenology has been explored for decades through intense experi- 
mental and theoretical effort. However, there is still little understanding of the connection to QCD and the basic build- 
ing blocks of nature, quarks and gluons. For instance, as a first "benchmarking" step, Lattice QCD should post-diet the 
large nucleon-nucleon scattering lengths and the existence of the deuteron, the simplest nucleus. The connection be- 
tween QCD and nuclear physics will be firmly established with Lattice QCD, and will allow for an exploration of how 
nuclei and nuclear interactions depend upon the fundamental parameters of nature. In particular, it is believed that an 
understanding of the fine-tunings that permeate nuclear physics will finally be translated into fine-tunings of the light- 
quark masses. While these issues are of great interest, and it is important to recover what is known experimentally to 
high precision, these goals are not the main objective of the Lattice QCD effort in nuclear physics. The primary reason 
for investing resources in this area is to be able to calculate physical quantities of importance that cannot be accessed 
experimentally, or which can be measured with only limited precision in the laboratory. Two important examples of 
how Lattice QCD calculations can impact nuclear physics are in the structure of nuclei and in the behavior of hadronic 
matter at densities beyond that of nuclear matter. Lattice QCD will be able to calculate the interactions of multiple 
nucleons, bound or unbound, in the same way that it can be used to determine the two-body scattering parameters. For 
instance, a calculation of the three-neutron interactions will be possible. 

Extracting hadronic interactions from Lattice QCD calculations is more complicated than the determination of the 
spectrum of stable particles. This is encapsulated in the Maiani-Testa theorem [2], which states that S-matrix elements 
cannot be extracted from infinite-volume Euclidean-space Green functions except at kinematic thresholds. This is 
clearly problematic from the nuclear physics perspective, as a main motivation for pursuing Lattice QCD is to be able 
to compute nuclear reactions involving multiple nucleons. Of course, it is clear from the statement of this theorem how 
it can be evaded, Euclidean-space correlation functions are calculated at finite volume to extract S-matrix elements, the 
formulation of which was known for decades in the context of non-relativistic quantum mechanics [3] and extended 
to quantum field theory by Luscher [4, 5], The energy of two particles in a finite volume depends in a calculable way 
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upon their elastic scattering amplitude and their masses for energies below the inelastic threshold. 



MESON-MESON INTERACTIONS 

Due to the chiral symmetry of QCD, %% scattering at low energies is the simplest and best-understood of the hadron- 
hadron scattering processes. The scattering lengths for nn scattering in the s-wave are uniquely predicted at LO in 
[6]: m % d^ = 0.1588 and m % d^ = —0.04537 when m K is set equal to the charged pion mass. While 
experiments do not directly provide stringent constraints on the scattering lengths, a determination of s-wave nn 
scattering lengths using the Roy equations has reached a remarkable level of precision [7, 8]: m^a 1 ^ — 0.220 ±0.005 
and mjfl'jj 2 = —0.0444 ±0.0010. At present, Lattice QCD has computed %% scattering only in the 7 = 2 channel with 
precision as the 1 = channel contains disconnected diagrams. The only existing nf = 2+1 Lattice QCD prediction 
of the n + n + scattering length involves a mixed-action Lattice QCD scheme of domain-wall valence quarks on a 
rooted staggered sea [9]. The scattering length was computed at pion masses, m n ~ 290 MeV, 350 MeV, 490 MeV 
and 590 MeV, and at a single lattice spacing, b ~ 0.125 fm and lattice size L ~ 2.5 fm [9]. The physical value of 
the scattering length was obtained using two-flavor mixed-action ^-PT which includes the effect of finite lattice- 
spacing artifacts to &{m%b 2 ) and G(b A ) [10]. The final result is: m % a l =£ = -0.04330 ± 0.00042. Lattice QCD 
calculations at one or more smaller lattice spacings are required to verify and further refine this calculation. Recently, 
the ETM collaboration has performed an/ =2 calculation of the n + n + scattering length at pion masses ranging from 
m K ~ 270 MeV to 485 MeV, at two lattice spacings and in two lattice volumes. The result extrapolated to the physical 
pion mass is: m^a 1 ^ = -0.04385 ± 0.00028 ± 0.00038. 

In Figure 1 (left panel) one sees that the n + n + scattering length is consistent with the current algebra result up 
to pion masses that are expected to be at the edge of the chiral regime in the two-flavor sector. While in the two 
flavor theory one expects fairly good convergence of the chiral expansion and, moreover, one expects that the effective 
expansion parameter is small in the channel with maximal isospin, the lattice calculations clearly imply a cancellation 
between chiral logs and counterterms. However, as one sees in Figure 1 (right panel), the same phenomenon occurs 
in K + K + where the chiral expansion is governed by the strange quark mass and is therefore expected to be much 
more slowly converging. The n + K + scattering length exhibits similar behavior when ytKn a K+n+ i s plotted against 
HKx/y/ficfn- This cancellation between chiral logs and counterterms for the meson-meson scattering lengths is quite 
mysterious. 
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FIGURE 1. 

wi7i®7i+ n+ wiji/fn (left panel) and iiij ( ci^+ vs. tn^/ ' fy_ (right panel). The black points are the results of Lattice 
QCD calculations by the NPLQCD collaboration, the blue point (left panel) is from the CP-PACS collaboration, the purple, red and 
magenta points (left panel) are from the ETM collaboration. The red lines are the current algebra predictions. 



BARYON-BARYON INTERACTIONS 

The NPLQCD collaboration performed the first n/ = 2+l QCD calculations of nucleon-nucleon interactions [11] and 
hyperon-nucleon [ 1 2] interactions at low-energies but with unphy sical pion masses, and the nucleon-nucleon scattering 
lengths were found to be of natural size. Such calculations are far more challenging than those involving mesons due 



to the exponential degradation of the signal-to-noise at large times. A summary of all Lattice QCD calculations of NN 
scattering is shown in Figure 2. The results calculated on the anisotropic clover gauge configurations are consistent 
with those obtained with mixed-action Lattice QCD [11], and it is interesting to note that the results of quenched 
calculations [13] yield scattering lengths that are consistent within uncertainties with the fully-dynamical iif = 2 + 1 
values. 
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FIGURE 2. A compilation of the scattering lengths for NN scattering in the 1 Sq (left panel) and 3 5j (right panel) calculated with 
Lattice QCD and with quenched Lattice QCD. The vertical lines correspond to the physical pion mass. 

Figure 3 shows a summary of the high-statistics results for baryon-baryon interactions recently obtained on the 
20 3 x 128 anisotropic clover gauge configurations with m n ~ 390 MeV. The AA channel is found to be negatively 
shifted in energy which may signal that the lowest-state is in fact a bound-state, but calculations in a larger volume 
are required before more definitive conclusions can be drawn. Calculations at this pion mass should be viewed to be 
only a proof of principle, and in order for lattice calculations to provide meaningful constraints on interactions at the 
physical pion mass through chiral extrapolation, the pion mass must be substantially reduced down toward the physical 
value (m K ~ 140 MeV) while maintaining the integrity of the calculation (i.e. small enough lattice spacings and large 
enough volumes). 
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FIGURE 3. Baryon-baryon interactions extracted from calculations at m n ~ 390 MeV [14]. 

I should also mention that there have been calculations of energy-dependent and source-dependent "potentials" that 
are constructed in such a way to reproduce the calculated phase-shifts at the energies of the states in the finite lattice 
volumes, for instance, see Ref. [16]. Currently unjustified assumptions are required in order for these objects to be 
compared with the modern nucleon-nucleon or hyperon-nucleon potentials that are used in the calculation of nuclear 
and hyper-nuclear structure and interactions, and I will discuss them no further. 



MULTI-BARYON INTERACTIONS 



The first significant steps towards the calculation of the properties of nuclei directly from QCD were taken by the 
NPLQCD collaboration [15] and the PACS-CS collaboration [17] during 2009. The NPLQCD collaboration performed 



2 + 1 calculation of the correlation function with the quantum numbers of the strangeness s = —4, baryon 
number B = 3 system, which is labeled as "nS°E " for convenience, and also of the correlation function with the 
quantum numbers of the triton (or 3 He), at a pion mass of m K ~ 390 MeV [15] as shown in fig. 4. Further, the PACS- 
CS collaboration performed a quenched calculation of the system with the quantum numbers of the triton and of the 
a-particle at a pion mass of m n ~ 800 MeV [17] as shown in fig. 4. 
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FIGURE 4. The left panel shows the effective mass plot in the triton channel resulting from the nt = 2 + 1 calculation with 
m K ~ 390 MeV after the subtraction of 3 x Mtf. The right panel shows the quenched results for the binding energies (in lattice 
units) obtained by the PACS-CS collaboration in the channels with the quantum numbers of the triton and the a-particle [17] with 
m % ~ 800 MeV. 

The dream of being able to perform reliable calculations of the interactions among multiple nucleons and hyperons, 
and of the structure and reactions of light-nuclei, directly from QCD is starting to be realized. The path forward is 
clear, and the next decade will be a truly remarkable period for nuclear physics. 
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to quantum field theory by Liischer [4, 5]. The energy of two particles in a finite volume depends in a calculable way 
upon their elastic scattering amplitude and their masses for energies below the inelastic threshold. 

MESON-MESON INTERACTIONS 

Due to the chiral symmetry of QCD, nn scattering at low energies is the simplest and best-understood of the hadron- 
hadron scattering processes. The scattering lengths for n% scattering in the s-wave are uniquely predicted at LO in 
%-PT [6]: = 0.1588 and mjs'jf = —0.04537 when m K is set equal to the charged pion mass. While 

experiments do not directly provide stringent constraints on the scattering lengths, a determination of s-wave nn 
scattering lengths using the Roy equations has reached a remarkable level of precision [7, 8]: m^a 1 ^ = 0.220 ±0.005 
and m^a 1 ^ = —0.0444 ±0.0010. At present, Lattice QCD has computed %% scattering only in the 7 = 2 channel with 
precision as the 1 = channel contains disconnected diagrams. The only existing nf = 2+1 Lattice QCD prediction 
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flavor theory one expects fairly good convergence of the chiral expansion and, moreover, one expects that the effective 
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FIGURE 1. m K a K + K+ vs. nix/ fx (left panel) and mi c a K + K + vs. m^/fy (right panel). The black points are the results of Lattice 
QCD calculations by the NPLQCD collaboration, the blue point (left panel) is from the CP-PACS collaboration, the purple, red and 
magenta points (left panel) are from the ETM collaboration. The red lines are the current algebra predictions. 
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FIGURE 2. A compilation of the scattering lengths for NN scattering in the 1 Sq (left panel) and 3 5j (right panel) calculated with 
Lattice QCD and with quenched Lattice QCD. The vertical lines correspond to the physical pion mass. 

Figure 3 shows a summary of the high-statistics results for baryon-baryon interactions recently obtained on the 
20 3 x 128 anisotropic clover gauge configurations with m n ~ 390 MeV. The AA channel is found to be negatively 
shifted in energy which may signal that the lowest-state is in fact a bound-state, but calculations in a larger volume 
are required before more definitive conclusions can be drawn. Calculations at this pion mass should be viewed to be 
only a proof of principle, and in order for lattice calculations to provide meaningful constraints on interactions at the 
physical pion mass through chiral extrapolation, the pion mass must be substantially reduced down toward the physical 
value (m K ~ 140 MeV) while maintaining the integrity of the calculation (i.e. small enough lattice spacings and large 
enough volumes). 
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FIGURE 4. The left panel shows the effective mass plot in the triton channel resulting from the nt = 2 + 1 calculation with 
m n ~ 390 MeV after the subtraction of 3 x Mtf. The right panel shows the quenched results for the binding energies (in lattice 
units) obtained by the PACS-CS collaboration in the channels with the quantum numbers of the triton and the a-particle [17] with 
m % ~ 800 MeV. 

The dream of being able to perform reliable calculations of the interactions among multiple nucleons and hyperons, 
and of the structure and reactions of light-nuclei, directly from QCD is starting to be realized. The path forward is 
clear, and the next decade will be a truly remarkable period for nuclear physics. 
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